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SYI’lOfBIS 


An application of FM technique to estiiaate the 
machinJjQg zom (or priiarj shear deforaation zone) dlaensions 
in orthogonal metal cutting has been cawied «it. In this 
study, ii»tal undergoing plastic defosMation is treated as a 
non-Bewtoniaii, incoH^ressihle fluid* Tbs analysis iirrolves the 
prediction of strain invariant ccmput'ed using ttie veloei% 
dlstrihution in Ihe cutting zom using fluid flow equations of 
mass halance and aoaentua halance* For the sal» of simplicity, 
it is assumed that Idxe cutting zone is isothermal and no 
tei^eretui^ effec'te on the properties of work material are 
included* Ihe non-lir^ar fluid flour gowinalig equations are 
solved iteratively wing the FM procedure based on the 
GalerklnO's approach, Ihe applicability of the model is 
confimed by c«paring the pres«it insults with the available 
experimental data from earlier stt^ies* It is eatablls,hed 
that a cut off value of 3 */• of tkm strain rate invariant 

can be used for approxlmtely demrcating the shear zone# 


CHAPTER -I 


1.1 GEMSPvAL BACKGROUICi : 

Metal cutting forms the basis of the modem 
engineering industry and nearly MO % o£ all the engineering 
products produced every year involve metal cutting directly 
or indirectly. A basic knOTledge of metal cutting process will 
be helpful in solving various problems encountered during 
saohinlng in a scientific way* Purther» a detailed study of 
the stresses acting during mchining could lead to the 
prediction of cutting forces aM ISie total pcwer input required 
for a certain rate of jKtal removal* The understanding thus 
gained can help in increasing the prcsductivily and lowering 
the machining cost anui also in tee development of new tool 
materials. Although a large voIum of tee experimental data 
about forces and stresses in i»tal cutting is available. There 
are not many teeor/tical models teich can ; accurately predict 
these variables in metal cutting# 

The stechanics of metel cutting is usmlly studied 

■.AJ'' 

in terms of tee formation of chip# The chip forma -Hon 
meteanlsm involves tee intense plastic deformatios of tee metal 

at stresses exceeding tee yield steess value at h%h strain rates, 

e 

It stron^y depeMs upon several pammewms such as tee cutting 
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o£ tile workpiece »aterial, Sssperiaental otjservations show that 

O' 

the shear zoae. In the chip fornation process, the friction on 
the tool chip interface also plays an important role* In the 


heen investigated in detail, Bie deformation process within 


1.2 RE¥iE« m mmious work s 


is tbs classical single shear plane model of Merchant and 
Piispanen According to this model, the formation of the 










evaluated the values of flow stross as a functim of strain 
fcB* a range of strain rates. 
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laWj Tifherein the flow stress is partitioned, into two parts, 
a non-theiml part which occurs in the shear fronts (or shear 
bands "i and a thewial nart which occnr' In the laaalla rations ^ 


He has also investigated the dislocation distribution %iftLlch either 
exists in the metal prior to cutting or is produced by Hie 
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(a) Solve Idae flow problem for a given temperature 
distribution , 

(b) Solve the thermal equation to obtain the temperatures, 
with a known flow field, 

(c) Iterate the floi#? solution with the plasticity vali^s 
adjusted according to the temperature field* 

At present, there are no FM studies of the nKstal- 
cutting process using -tiie viscoplastic flow approach, 

1*3 OBiTECIIVES AHD SCOPS OF iHE PRESENT STUDY i 

In the present study, the visco -plastic flow approach 
proposed by the earlier authors is employed for 

investigating the metal cutting problem* For the sake of 

i simplicity, the thermal effects (which, are in fact quite 

' ' < 

: important for evaluating the material properties) are not ^ f 
considered and properties are evaluated at a topical averse 
temperature of the cutting zone, Ihe plastic flow of the metal 
is modelled as a highly viscous, non-Uewtonian flow with viscosi 
values strongly depending lapon the stiraiin rates* Considering ’Sue 
plastic deformtion processes to be doisinant in the shear son® 
and -assuming incouqpressible, pure shear deformtioa, the resultl 
flow equations are solved with the help of the finite Element 
procedure • 



23ae objectives of the present work are to determine 
the strain rates in the cutting zone and to establish a 
theorltical model for the plastic deformation processes which 
eoB^sares favourably with the known experimental results* It is 
also desired to predict the shape and size of shear defoimiation 
zones, based on suitable criteria with regard to the 
bouMarles . 
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CHAPm-II 


FLOW MODEL FOR METAL CUTTIHG 


A two dimensional orthogonal machining uncer dry 
cutting conditions in considered as shovm in Pig, (2.1), Ihe 
tool is assumed to be perfectly sharp with no built up edge 
(SUE)* She work material is taken to be ductile resulting in Ihe 
formation of a continuous chip, Ihe tool material on the 
other hand is taken to be perfectly rigid without any deforma- 
tion* In the present work, metal cutting is assumed to be 
steady and isothermal. A flow stress model based on the 
plastic behaviour of the work material is developed here to 
investigate the various features of the metal cutting process, 

2,1 Ifechanica of cutting and chip formation j 

Ihfi TOtal cutting process can be analysed by s tucty. ing 
Idle formation of chip in detail. As uncut metal approaches laae 
tool, it is severely compressed in a region in front of 1die 
cutting tool, %ils causes plastic flew of the metal accompartaed 
by high temperature if the metal is ductile. When the stress in 
wcark^^iece ahead of the cutting' tool reaches a value exceeding 
the ultimate strength of the metal, the particles shear to form 
the chip which moves up along the rake face of the tool. She 
outward or shearing movement of each successive element is 
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Tool 


c=> 

Workpiece 


Fig. 2.1. Geometry of Chip Formation in Orthogonal 
Cutting. 
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arrasted by work hardening and the moveiaent is transferred 
to the next element* Ihis process repeats and a continuous type 
of chip is formed j having a highly compress oa and bu.r:*5 
underside and minutely serrated top side caused by the shearing 
action* 'iiie plane along which the elements shear is called 
the shear plane. 

Ihe deformation does not occur sharply across ^ the 
shear plane j rather it occurs in a narrow band, Ifetallographic 
studies (5) have shown that the grains of chip material begin 
to elongate along the line AB below the shear plane .and 
continue to do so until they are completely deformed along ihe 
line CD, above the shear plane as shown in Fig, (2,2). fhe 
region between the lower surface AB, and the upper surface CD 
is called the shear zone. It has been noted 
shear zone itself consists of two subregions, one in ^aich the 
metal shears parallel to the shear plane and another in which 
shearing takes place parallel to die tool surface, 

fhus the zone in -s^&ich plastic flow takes' place parallel 
to the shear zone is known as the primary shear deformation 
zone (PSDE) . %a primary zone has been observed to be wedge- 
shaped, and its thickness depends upon machining parameters 
such as cutting velocity. For , instance, it has been reported 
that an increase in the cutting velocilgr causes -the shear zone 
thickness to z*educ# Plas'tlc deformation of the metal 

occurs not only In Ihe prima^ry zemet but also In a region close 
to tool-chip interference. As «i« Chip moves over ihe 
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too^ace^ a very thin layer of chip iisSiich is in contact with 
the rake face of the tool undergoes additloneil plastic 
deformation* Ihis phenomenon is called as secoMaiy shearing 
and the zone in ^ich secondary shearing occurs is known 
as Secondary Shear Deformation (SSDZ) Zone . In the secondary 
Zone, each successive layer of metal moves at a different 
velocity because of the plastic flow, fhe machined surface 
below the flank face of tool also undergoes a certain 
deformation* fhls happens because of the elastic after effects 
of the machined surface and flank wear of the tool lead to 
the high friction between -the mchined surface and tool flank, 
Thtis, modelling of the chip formation process requires a 
thorough imderstanding of the various plastic and elastic 
defonaation of the work material around iiie cutting edge, 

2,2 Iheory of Plastio Flew j 

Metals deform elastically only in a certain range of 
small strain* When a metal is strained beyond Its elastic limiti 
it starts yielding and some permanent (plastic) deformation 
takes place* If the metal is ductile, it undergoes a considerable 
amount of plastic deformation due to yielding until fj^cture 
occurs eventually* fhe behaviour of plastically deforming 
materials is rather complicated and fully rigorous* Accurate 
theories of such behavimir of arterial are not yet available* 
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Among plustic solids, an ideal plastic obeying 1iie Yon-mser-. 
yielding criterion gives a reasonable abstraction of the 
behaviour of certain real metals yielding plastically* Herei 
we shall briefly describe the basic constituents of the 
plastic flow theory of Von Mises and adopt it for the modelling 
of the metal cutting problem. 

Let us define the stress deviation 0^^ and the strain 


deviations ej|^^ at a point as ; 


°ii “ °Li'h °ii 

(2.1 (a) 

®id “ - 7 

(2.1 .a) 




Vhere 044, e » 4 are respectively the components of the 
stress and strain tensors and •y are respectively 

the mean values of the normal stress and normal strain* If 
the material is isotropic and obeys Hookers law upto the 
elastic limits then, 




( 2 . 2 ) 


where Q Is the modulus of rigidity. Beycnad Ihe elastic limit, 
e. . is no longer given by the Hookers law* In ihds case„ w® 
define the plastic strain deviation m the actual strain 
ffl’timigs the st3ra.in deviation that would be computed If the Koote’^s 
law is still applied. Let the plastic strain deviation be 
denoted by 



If at each instant of time, tlie rate of plastic 
deformation is then the increment of plastic strain in 

the time-interval dt is aM liie total plastic deviation 

after successful stages of yielding will be the algebeVaic 
sum of deformation that occurs at all stages, 

‘‘•U J 0 

where is laie initigl value of a|j^at time twO 

The following assumptioias are invoked in the modelling 
of a ¥on-Mises type ideal plastic solid* 

(i) Hooke's law holds for the mean stress and mean strain 
at all times. 


"u ■ ‘ii 


(2*6) 


Hence, the plastic strain is incompressible 


! as plastic 
Jn places where 
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(ii) Bie material Is elastic and obeys Hooker's as long as 

2 . 

the secoM in variant less than a constant K * 

Inother words ji no plastic deformation can occur as long asi 

jg *a ^ (2,7»a) 

and 0 when J^< (2,7,b) 

ij 

(iii) Yielding can occur (elastic limit is reached) only 

o o 

when J 2 » K . 'feJhen the yielding condition JTg-iC «0 
prevails, the rate of change of plastic strain is 
proportional to the stress deviation 

h “id 

Where, f* is positive factor of proportionality -^rtiich 
has the dimensions of coefficient of viscosity of fluid, 

(iv) Any stress state corresponding to ^ 
realized in the material, 

Ihe set of assumptions mentioned above contain two 
essential parts, namely the criteria for yielding and the stress - 
strain relationship in the elastic and plastic regions, Ihe 
yielding ©'■••dition is based on the second invariant of the stress 
deviation tansor. Such a yielding criterian was first proposed 
by Ym’-Mism ♦ aiie constant K can be identified with the yleM 
stress in simple shear. 
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We notice that since only the rate of plastic strain 
is specified hy the ah eve equa- ions fesr a varying loading 
sequence y successive plastic increments must he added together 
algeherically , Such a theory is called an Increitental Theory. 
Since the material behaves like a viscous fluid \fdien it is in 
the plastic state. It is called an ideal viscoplastic material. 

In many applications of theory of plasticity, the 
rate at ijdiidh a plastic flow occurs is of little interest. 
Sometimes an engineer is concerned only with the total amount 
of plastic flow. For such application, v/e often replace 
equation by the increment law - 

e^P) m K 0j. (2.9) 

ij 

Where k is a constant factor of proportionality and 
not a characteristic constant of matejfiflal. The sign of X Is 
detenained by the fact that plastic flow involves dissipation 
of energy iinder a given set of loading conditions, the value 
of X and hence the total plastic flow is determined by the total 
work done by the external flow* 

2,5 Behaviour of Hlastlc Material i 

In large deformation process such s’s machlnlEg involving 
plastic or Viscoplastic flow, the elastic deformation can be 
considered as negligible, such problems^ we find that the 
ccaistitative relation can be espresso in an Eulerian X&tm 
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linking the stresses and current strain rates# 'Uie situation 
is thus identical with iiiat of a flow of -riscouSj^ non-lfewtonian^ 
incosipressihle fluid Fig# 2,3 gives a partial classification 
of substances based on their rheological (i*e* shear stress Vs 
rate of strain) behaviour* 

r^s seen from the fig* non-Hewtonian fluids do not follow 
Ifewton®s law of viscosity (3hear stz^ss « cons*cant x shear 
strain rate) aM the assumption of a constant viscosity 
independent of teE^jerature, density ani shear strain rate does 
not hold good. 

Certain fluids, known as Bingham plastic fluids, 
behave as a rigid solid until a certain level of shear stress 
is attained anl behave like Newtonian fluids afterwards, as 
shown in Fig, /•2-3J7. A general plastic material will have 
a non-linear variation of shear stress with shear strain rate 
beyond a critical stress value (yield stress). 

2,4 Hew Stress Modelling % 

In this study > we shall model Itie metal undergoing 
plastic deformation as a non-Newtonian, incompressible fluid 
(no density change)' on the basis of the discussions presented 
in articles 2,2 and 2*3, 

For such a viscous and iacoippsssible fluid, we 
can write the constitutive relationship in the form 

^i3* ** ^iiJ^ ^ii^* i**4 (2*10) 

*•0 i 


where (-p) is -tixe mean normal stress and is 
the strain rates. 

The strain rate components in two dimensions can he 
written (Eqn,(2,11)) in terms of only the first (linear) 
order components of the deformation rate tensor^ neglecting the 
higher order velocity derivatives* 

®11“ li ' ®12 ®21 I ^ H ®22 “ ly (2.11) 


vhere u and v are the velocity components in x and y directions » 
respectively. Rewrite the strain rate tensor as s 


1 1 a 

w 


( 2 . 12 ) 


in which d^^ are the deviatoric components of stress , 

In the present work, for a Von-4tises type of plastic 
material modelled as an isotropic, in compressihle, non-Newtonian 
fluid, the variation in viscosity Is given hy (14) 

/ ^ Ar nt 

a * {.eh fs ) 

JX a " (2,13) 


/3 8 


In the above eqmtion, ju is the non-linear viscosity, 7 
and n are the physical ocnstants determining the vlseo -plastic 

heteviour of material aM Oy 1, th, antottal ylaXd .tr«,. 
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of tlie mterlal v^^ich is in general a function of temperature# 

* 

pressure and deviatric strain invariant e, fhe strain in 
variant can be represented in terms of the velocity derivatives 
as follows i 


& 

e 
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*2 

+ ®22 "• 


2e, 


12 ®21 



(2*14) 


2,5 Governing Sgtuatlons s 

For steady orthogonal metal cutting operation with 
the width of cut much larger than l^e feed or depth of cut# 
a two dimensional model may be employed, Treating the model 
as a non-Hewtonian# incompressible fluid iii the maidalning zone# 
the governing equations are given by i 

7 .V * 0, (mess balance for incois^ressible fluid) 

ea* •*> 

(2.15) 

Y* 7 ¥ « V. g (M(mentum balance) (2.16) 

In the above equation stress tensor can be written as s 

.7 

I #» -.p «t -P J +d (7 Y + T Y*^)- 2/5 M i 

(2,17) 



1 p 

J 


'I 



where p is the press\jre and d is the deviatorie stress tensor. 
Using equation (2,17) in equation (2,16) 

PI-SI = - V 2 I. (» {V I* I ) (2.18) 


In component £orm> from eqn,(2,18) we can writ# x and 
y momentum equation as 
X-momentura 


P Is + * ly I** <ly H >1 

(2,19.a) 

y-momentum 

P li - ^ “If * fe 

(2.19.b) 


I3ie mass balance equation (2*15) can be rewritten as t 



(2.19.C) 


2,6 BouMarv Conditions t 

I5i© boiiMary conditions for solving equations 
(2,l9#a,b,c) are basically of two types. Or siirfaoes 1,11^111,, 
? and ?I (Fig# 2.4) the velocity comp{»«tB are specified# 

On surfaces 17 aiad 7II aero traetloa aM «ero aoi»al velocity 
S|»oified. In mathemtical form# w« bave 
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(i) 

(ii) 


u V (cutting velocity) 

V a 0 on surfaces Ip II, III 

u a ¥ sina 

on surface ¥I 

V » V„ cosa 

c 


(2.20 .a) 
(2*20 .b) 


v/hers is chip velocity and a is the «ngle , 

V 

(iii) ?or surfaces r/ and ¥11. 


a 0 

n 


n ,0 at 0 

S3 S3 


(2.20.C) 


where v„ and n are normal velocity and the unit normal vector 
respectively at a point on the surface. 

(iv) On surface ¥, upto a length h, the slip velocity at the 
interface between the chip and the tool is presented* 
From the ejsperimental measurements of Tay et *al 

^s “ 5 ^c /“TT'WTir' 


V&ere ¥_ is the chip velocity and s is the distance from 
the cutting edge ireasured along the interface, ^yoM the plasti 
contact length h, the velocity is taken to be equal to the 
chip velocity on the surface ¥. Therefore* 

u « ¥^ J'Tlis/h sina 

A I for s<rh (2»20*<i) 

V « V^ J 1 +Bs/h COSEt 
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Fig. 2.3. Rhelogical Classification of Flow. 



Fig. 2.4. Problem Regions, Showing the Boundry 
Conditions. 
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® sina 
V =« Gosa 


for s:>-h 


(2*20.e) 


2,7 Hon~£)lmensloQalizatlon j 

I’/riting the Eons, £'ZA3{Q.f'b^Ci)J in a dimensiomless 
form facilitates generalisation to model a large range of 
problem. In Hiis respect, following dimensionless variables 
can be used 


s# x/t » y* = y/t 


u « u/y 


V* « vAr 


where V is the cutting velocity and t is the thickness of 
imcut chip. 

How, the dimensionless form of the continuity equation will be 


'du 


av^ 


m 0 


(2.21 ,a) 


writing liie dimensionless form of the x-momentum of equation. 



Let P*. p/Oy , n- p/j ^ 



vAiere is the uniaxial yield stress. 
On simplifying. 


0 3y 


* 


ax 


ax 


(2u 


* 3u 




dx 


''By 3x 

Siniilarly the ynmomentum Equation can be put in dimensionless 
form as i 


h (u"" K K) 

G 3x ay 


|E -j. (M^ + 

y ax ^ ay ax 


)1 


ay 


■* 

(2h (2*21 .c) 

ay 

From Eqn, (2.13)4 %e noimjalized. Tariation of 
Viscosity is given by 

, 1+ 

h » u/byt/V » (2*22) 

where B* » |p* ( 


rrtT 


The boundary condition can also be easily non- 
diiiensionalized as shown above# All velocities in Fig* (2*4) 
are scaled by ¥ and lengths are scaled by t* to obtain the 
dlaeasionless boundary conditions* 
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CHAPTER^ III 

FINITE 5LEM5NT AIIALYSIS 

3.1 INTRODUCTION t 

The solution procedure that has heen adopted for 
the prohlem under study is the Finite Element Method. Some of 
most Important advantages of the FEM procedure are the easy 
handling of a complex geomatrical shape of timt solution domain 
and the general manner in which boundary conditions can be 
introduced, father than requiring every trial solution to 
satisfy the boundary condi-*;lons, the condition may be presented 
after obtaining the eilgcberaic equations for assemblage. Sine# 
the boundary conditions do not enter into -the equations for 
the individual finite elements, one can use the same field 

variable model for both internal and boundary elements, 

0 . 

Morover, ihe field variable models need not be changed, when 
the boxuidary conditions change, la addition to the easy 
accMiodation of complex geq;f*ttry and boundary coxaHtions, the 
iEM technique has also been successful in representing 
co:^licated material property variation tl^t are difficult to 
inceo^car'ate into other numerical melhods for example, 

e 

formulations in solid mechanics 3have been devised for aniatroplo, 
n<na-linear, time dependent or temperature dependent material 
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i 

! 

behaviour* In "the present problem, the geoiaetry Is coiapllcated 
aai a strong non-linearity is present due to the high variation 
of viscosity in the plastic flov/ region. Iherefore^ the FE!>I ; 

procedure is quite suitable for obtaining the velocity field 
solutions, 

3,2 DSSCRIPTIOH OF FSM f2CHKiaUE t 

In PEM, the governing differential equations are 
converted into an equivalent set of integeral equations which 
are minimized in some sen, re over the whole solution domain. 

%ie solution domain is divided into several small elements of 
known shape ani the solution variable is approximated within 
each of these elements by suitable approximating functions. The 
integv^ral governing equations are then evaluated within each 
elemeiit using the approximating functions and all such elemental 
inte: .^rals are assembled to obtain the final matrix equation 
to be solved. The above steps can be summed lap as follows i 

(a) Obtain the equivalent Inte ,;ral equations for describli^g 
the problem. The boundary conditions of the problem 
would be transformed into appropriate boundary integrrals, 

(b) Divide the solution dcmaain into several element' and 
place chosen number of nodes within each element. Obtain 
the approximating functions for describing the field 
variable, within each tlwent in terns of the nodal valws. 
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L (/) = R (5.2) 

..'here, R is the residue which may be a function of the 
spatial coordinates and time* In Galerkin*s procedure, the 
residue, weighted appropriately by suitable weighting functions, 
is minimized over the whole solution domain D in an intc ral 
sense. In mathemcitical form* this becomes ; 

C Vk R dvw f Wj L (0*)dv » 0, for each i (3.3) 

where, Wj, are the weighting functions, 'Jhe approximation 0^ 
may be written as a piece continuous profile for each element asi 

i«1 ^ ^ 

In the above equation, » is the number of nodes per 
element and N^^^are the approximating functions for representing 
0* within the element (e) vdilch are called as the Shape Functions 
of th® element (e). TbB nodal values of ^ tor the element (*) 
are denoted by In Gal8rkin*s procedure,, the wei^ting 
function are taken to be the same as the shape functions 
within each element* 

3»3 nmB OF FORMULAflOH t 

fhere are four tj^pes of formulation i^feich are in use 
for obtaining the finite solution of any flav problem, fhese are t 

(i) sOTm mmnm 

(ii) im&M mmtim ahd ramuuxKM 

(iii) va^ocixr fressqbi raisjuxicif 
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(Iv) PEMLITT FUKCTIQE PORMULATIOH 

!Ehe stream function approach for solving a flow 
problem was first proposed by OSLON By introducing the 

stream function f, the continuity equation can be satisfied 
exactly and the momentum equation can be combined into a single 
function in terms of 'fhe main disadvantage of this formiulation 
is that -the governing equation is a fourth order differential 
equation and interpolation functions of high order are required,' 
Also in order to obtain velocity and pressure, further 
processing of results is required* 

In the stream function vorticity formulation, in stead 
of dealing with a single governing equation fourth order, the 
stream function (f) and vorticity (w) are chosen as the unknowns 
and two governing equations of second order in and w result 
riuT. In this formulation, the boundary couiitions may consist 
of specified values of f and w or specified values of their 
first derivatives. Ihe principal difficulty in using the stream 
fimcticai-vorticity foimiulation is that, in general, vorticlV 
is unknown a priori along solid boundaries. 

In the velocity pressure formulation, the pressure and 
velocity components are considered as the nodal variables, 

Ihis formulation isused because of the following advantages t 

(1) the formulation is readily extended to thir^ee dlaensioiw# 

(2) toy C° (linear element) continuity is required of ihe 

element interpolatltm Idncttoas* 



(3) Pressure, velocity, velocity gradient and stress 
boundary conditions can be directly incorporated into 
the matrix equations : 

(4) Rree surface problems are tractable . 

(5) The formulation requires less computational time than 
stream function and vorticity fonaulation (17). 

1 

The penality function approach ifleJZ is #ery populair 
procedure in viiich the pressure variable is eliminated and only 
the velocity components are solved* In the present study> we 
shall use this approach because retaining pressure as a variable 
often leads to numerical difficulties. The penality formulation 
involves treating continuity equation as a constraint among the 
velocity components and the continuity constraints is introduced 
into the momentum equation in a Lagrange multiplier K, It is 
possible to choose \ to be very large and eliminate the pressure 
in terms of X. Ihe boundary condition to be specified for this 
formulation are either specified values of velocity components 
or specified values of their first derivatives. 

3.4 DERIVATION OF FINITE ELEMMT EQUATIOHS : 

Ohe dimensionless governing equations (after dropping 
the stars for convenience ) are it 


I (« || + v||) * "H + |l (2fx ||) •*- If)} 
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( 3 , 5 * 13 ) 


3u 

3x 



(3,5 ,c) 


In the present wcrk^ the inconprassibility constraint 
i:Iqn,(3,5,c) has been treated with the help of a penality function 
approach. Introducing a parameter K such that X>>pj the 
in compressibility condition can be approximately written as t 

H + ^ - pA (3.6) 

Wo¥j using Eqn,(3,6) in Eqns.(3,6,a and b), we obtain : 





( 3 , 7 ,a) 


1 


0 





(3,7,b,) 

The unknown u and t in the aboye equations can now b« 
expanded in terms of ncdal yalues and appropriate interpolation 
function in the form i 





u 


a 
a S 
i«1 


*’1 


% 


a 


IT « S 




(3*8,a) 

(3.8,1>) 


wherei and are the interpolatir^j functions (which need 
not necessarily be of the saras order) and Uj^ and are the 
nodal values of the unkno’.ai* S’ollowin" Tanada et# al 
in order to simplify the algeberaic structure of the equations^ 
we select « mJ# Ei^t noded isoparametric elements are used 
for dividing the solution domain, For the isoparametric elements, 
the global coordinates x,y of apoint v/ithin are element may 
be transformed as j 
8 

X w s Xi (3.9,a) 

i i 

8 

y = E N, y. (3,9,b) 

is»i 

fhe velocity shape functions, for the for 6 -noded 
isoparametric element are given in Fig, (3,1) in tsrro of 
local coordinate 5 , ^ by the following foCTaulae s 
Corner ncKies *, 

% « I (1+% S) (% h~1) 


« 


(3.10 ,a) 
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Mldside nodes : 

N » 1 ( (5*10. b) 

««« £t* JL 

% “ ^ (1+^i i) (1-^^) (5.10.C) 


Employing the Galerkin*s weighted residual approachj 
Equation (3»7,a) becomes. 


ne ^ ^ /n n 3H^ n 

2 j I lx ^ 

K 1 y + 2 IHA ^ ^ 5 f2u 2 

^ I ”3X V i 3y 1. H ^ 

rn s ^ u^4. r ^ AA « 0 



(3.11) 


Vi/here the outer summation is overall the area ele»eats 
in the domain end the inner siimmation is overall the appropriate 
number of nodes of an element. In the present ease, the number 
of nodes per element, n»8# 


Invoking the Green's theorem , the second order texw 
can be partially integerated as * 




aM 




, e 

ii 


fy0 NlZ~<2ll ||) * H) "yJTil! 


f ^ 

^ ai‘L n aJL 

rr •^"- 1 “Sx "a 


•5* 




’.’i/here denotes the element surface, \idien the 
contribution from two adjacent elements are sirmaed up a shown 
in Fig* (3,2)* fhe find contribution from the side becomes zero 
unless it happens to be part of a boundary, 

Tlie required form of Equation 3.11 is 

I o(% \ “fc **k * ‘Sx "a 

8H, 3H. 3N, 8H. 3N. BN 3H, 3N. ^ ^ . 

* W‘ 8y V "W “ « « “a* ^ 'Sy ■»’a)-/“ 

■ 4® % • 0 

(3.12) 

In the above equation, the surasation sign I has been droppf 
by adopting the convention that summation is ii^lied if a 
subscript is repeated In the saise form. 
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KoWj tlie Qomentura equation in y direction may similarly be 
v/ritten as s 




(2?i|^) Hy+X (|H |X) 


(3.13) 


Hie final assembled form of the aboye equations 
(3,12# 3,13; i& 

where, the vector x c entrains all the nodal unknown u^ and v. , 

It is possible to write the coefficient matrix ^ 

sub^trioes r= iiJ. as follws : 


^11 ^12 

Cgi G22 J 


dA 


ne /" 

/"A 7 * S 1 
i, I j^e 


(3.15) 
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where 



In the above expressions^, we Mv@ treated u and v as 
Imown values (froa previous iteration) for the sake of Mklng 
the matrix equation (3 *14) qmsi*linear* 


On the boundaries where natural bouaiary conditions 
are imposed^ these line integerals appear in the ri^t haM 
side vector as illustrated belcM : Ih© oolxima vector is 
written as i 

hi 

\ da (5.17) 

bgj 


ne 

S 

1 


/. 


•fehere 

b^ -» £*(2h *** H) ^ -7 




C3*l8*a) 


= h ^ i )«, -^(2^ If) (|S ^ f ) “y J 


(3.l8.b) 


The matrix equation (3,14) is non-linear due to the 
variation of viscosity p and also due to ineratial terms. 

Therefore, one needs to employ a suitable iterative procedure 
for determining the nodal imknowns. 

3*5 mTRIX SOLUTION TECHMIGUS t 

The solution procedure used her® for the solution of 
the assembled matrix equation is known as -tiie Frontal Method 
Z"20_7‘. When tbe number of unknown in the problem is very large ^ 
the method adopted fur solving the assembled equation has a 
significant bearing on the computer storage requirement and 
execution time. It is very difficult to solve these equations 
by simple matrix inversion techniqms. The Frontal Method is 
based on direct Gaussian ellimination for solving the 
symmetric matrices where the leading diagonal Is always used 
as apivot, For unsymmetric metrices, encountered in a wide 
range of engineering problem, the most suitable pivot is not 
necessarily on ttie leading diagonal and frontal solution routine 
exist for off diagonal pivotiiog (6), This however tend to be 
more time consuming# The method used here uses only diagojaal 
pivoting incorporates many features of the Frontal' method for soXvii 
sywetrio matrices* 
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In general terms > the overall solution technique 
consists of 5 

(1 ) Formation of Element matrices 

(2) Assembling into a Global matrix 

(5) Introduction of BouMary Conditions 

(4) Reduction of Global fluid matrix using a Gaiissian 
elimination technique, 

(5) A back substitution process, 

3,6 BASIC PHILQSPHY OF TBB FRQHIAL ICBTHQD i 

Ibe primary objective of Idie frontal meliiod is the 
elimination of variables as soon as possible after their 
introduction, into the global matrix via the appropriate 
equations > As soon as the contribution from all the element# to 
a particular nodal point have been assembled, 1iie corresponding 
variable associated with liiat node can be eliminated, lEhe 
complete matrix is therefore never assembled since the reduced 
equations can be eliminated from cor® and stored on disc# Ibe 
equations held in core, with the corresponding nodes and 
variables are collectively called as the Front and the number 
of unknown variables in the front is termed as the Front 
Ihe front width changes ccmtinually since only after the 
contribution to a particular node have been fully summedt the 
corresponding equation reduction based on a diagonal pivot 
can be executed. 


A preassigned global matrix core area is first filled 
from contributing elements and then the largest diagonal entry 
in the preassigned core is foimd and used as a pivot in a 
direct Gaussian Elimination process « As the maximum pre -determined 
number of equations are eliminated the corresponding reduced 
equations are written on to the disc and moi^ elements and 
corresponding equations are taken into core. ?he Fig, 
shows the basic idea of the Frontal Ilethod, The equations^ 

nodes and variables currently core are termed Active j, 
those assigned to disc are D eotivotea and those yet to 
appear in core are Inactive . 



Fig. 3.1. Eight Noded Element Showing the Local 
Co-ordinate System. 



Fig. 3.2. Typical Element Interface Where Residual 
of Surface Integral 0. 



O Inactive Nodes 
• Active Nodes 
0 Deactivated Nodes 


V//X Front 

^✓✓^J ^Assembted Element 

I— — n4 

I • ‘ Next Element for Assembly 
Fig. 3.3. Dafinitions of Front and Elamant Numbering 
for Minimum Frolit Width. 
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CHAPTER «iy 
RESULTS AHD DISCUSSIOH 

Using WWl technique in vlscoplaatic model of metal 
cutting^ results are obtained for strain rate invariant, stream 
lines and mean width of primary shear deformation zone. To test 
the validity and accuracy of the mchelg analytical results 
are compared with essperimental results. To test the sensitivity 
of the model, a parametric study Ihs also been conducted. 

Results and discussion have been divided into three subheadings, 

(i) Mean width of primary shear deformation zone 

(ii) Flow lines during metal cutting 

(lli)(a) Effect of different parametr^s on the mean width of PSDZ 
an- As/t . 

(b) Dls ia^ibutlon of strain rate invariant In PSDZ. 

(1) Mean width of Primary shear deformation zone t 

In Fig, 4 ,1*4,3, constant strain rate invariant 
curves have been plotted for three different cutting conditions . 
a© strain rate invariant (5)has been noraalised with respect 
to its maximum value ( each cutting condition, 

graphs of normalized isostrain rate ^variant (e/Saiax^ 
to 3 y * 4 */» 5 V* have been plotted* 
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In order to estimate the size rSDZ (or mean width 
of PSDZ)s normals exterdlng up to the desired mrve say 5*/, 
cutt off pointy are drawn on liie shear plane (Fig, 4.1 )« Ihen 
the average of all these normals is calculated and is defined 
as the mean width of PSDZ It is to he noted that secondary 

shear deformation zone and a part of PSDZ ( shaded area for 5 */. 
curve in Fig, 4,1) is excluded to simplify the measurement 
procedure. However, better method to determine As, could be 
to divide the total area imder the curve by the shear plane 
length. After comparing the mean width of PSDZ with the e^tperiaental 
results (Table 4*1), it Is seen that for all he three cutting 
conditions, the three percent curve consis ..ntly gives the 
primary shear zone thickness '^ich agrees will with the 
experimental results, %erefore, it was decided to use the cut-off 
value of normalized strain rate invariant as three percent in 
the present analysis for investigating the properties of th# 
shear zone* 

33ie shape of thre^ percent curve iraiicates that the 
primary shear zone is approximately Wedge-Shaped and the 
secoi»iary shear zone ad;5acent to the rake face is approximately 

triangular (shadow portion near the tool.ohlp ooniaot). Also 
It IS seen (Fig. 4.1-4.3) ttiat the plastic deformation aone 
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TABLE 4.1 


Comparision of Experimenlial ( Ref, 6 ) and Theoritlcal 
value of Biean width of Rrisiary Shear Defonaation zcmm. 


S*No 

CuttiUK 

Condition 


Mean width 

of PSDZ ( 

s/t) 


(Degree) Degree 

u 

(m/s) 

t(mm) 

E:^eriaental 

cut-off value 1 







4/, 5V* 

1 

10 

25.0 

5.07 

1*75 

0*725 

0.743 

0.686 0.410 

2 

10 

28*0 

2,05 

1 .05 

1.05 

1.0S3 

0,820 0.624 

3 

10 

28*0 

4.15 

0,648 

0,648 

0 ,668 

0,459 0,382 


{ 

i 



In Fig* the flow paths of metal particles 

during metal cutting are presented for two different cutting 
conditions. It is seen f3s»om these plots that the paths of the 
particles are straight lines in both ‘the uncut metal and fully 
f canned chip* and sharp turning takes place in a thin zone in 
front of the cutting edge. It is this turning that contrihxrfees 
chiefly to the shear strain rate and the plastic flow of the 
metal occurs in this zone only, jPhotographical observatioas of 
printed-grids over the workpiece by Stevenson and 0:^#y 
have shown that the flow path of metal particles can be 
approximated by hyperbolic strain lines, Ihe results of the 
present study are in good agreement with these experiiwntal 
observations* 

(iU)(a) Effect of different parametres o n the mean 
PSDZ and AS /t t 

Having established the validity of 'te.eoritieal model* 
a set of hypo-iittitBa cutting cosditl« 0 » were considered in 
order to study the effects of iSm i^rtant cutting paremetres 
such as cutting velocity , (v)# (t) toe toear angle 
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upon the defopmation zone thickness, !Qie effect of the cutting 
velocity (v) upon the dimensioi^ess shear zone thickMSS ( 
in shown in Fig, 4,5, It is seen that ^^a/t decreases with 
Ihis trend has also been observed eagperimentally by Kakayoma 
and Okushiiaa ani Hitonii • From the theoritical point of vie^t 
the decrease of ^s/t wilii V may be explainea as follows. It® 
extent of the shear deformation zone depends upon the value 
viscosity. For a low value of viscosity * a smaller zone is 
and vice-versa, A higher cutting velocity (at const t) results iu 
a higher strain rate, which leads to a lower value of viscosity us 
seen from equation 2,13, therefore, at higher speeds > the 
dimensionless shear zone thickness is smaller. 

In fig, 4,6, the effect of the shear angle upon 
is shown. It is seen that Zss/t decreases with 0, For small®** value 
of 0, chip becomes thicker and so the length of the shear ^ 

deformation zone Increases* However, if the amount of matei*!^ ^ 
be removed and the rate of metal removal are fixed (i,e, 
value of t and v), the total deformation undergom by the mel^ 
is more or less constant* Therefore, while the length of the shear 
zone increases, due to an increase in 0, its thickness reduce* • 

The trend seen in the present result was earlier noted by Fan<i«y 
and Jain (23), 

The effect of the undeformtd chip thickness t, upon 
mean width of Ihe primary shear zone as is shown in Fig* 

The mean width A s increases approatimtsay in proportion 
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as was noted by Nakaynama _£*22j7. It appears that angular sha|« 
of 1216 shear deformation zom is more or less a constant when all 
the parametres other than t are constant* therefore* as increases 
with t. 

(iii)(b) Distribution of Strain Invariant in PSDZ t 

In Pig* 4.8, the variation of the a dimensionless 
parameter £H* Cl/(v/t) ) J is plotted against the oiher 
dimensionless parameter /“distance i /tJ7 measured along the shear 
plane from the cutting edge* As expected, high strain rate values 
occur in 1h@ vicinity of the cutting edge and far away from the 
cutting edge, f* is smaller* In metal cutting theory, for the sake 
of simplicityi it is often assumed that the shear stress and r&im 
of the shear strain are constant along the shear plarie* Ihe present 
results show that the strain rate is not constant and that a 
better assumption be that shear strain varies with i (Pig* 4*8)* 

In Pig* 4.9, 12ae dimenalGnlssB parameter is plotted 
against the dimensionless distance (h/t) travelled by a partldle 
through the deformation zone. It is found -that ^ decreases on 
either side of the shear plane as the particle moves away from -tee 
shear plane * It is obvious because on gw side It is moving towards 
the uncut work material where e* is apprcaiiwtely »aros On the 
other side, it is moving towards ’fee felly forwd chip where no 
further deformation takes place, h«n»ti ^ wd^d again |ttftin a 


zero value. 




4 % 

5 Vo 






Fig. 4.1. Constant Strain Rats Curves During Metal Cutting. 
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Fig. 42. Constant Strain Rats Curvss During Mttai Cutting, 
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Fig. 4.3. Constant Strain Rate Curves During Metal Cutting. 



Fig. 4.4q. Flow Lints During Metal Cutting. 
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Fig. 4.6, Effect of Shear Angle on Dimensionless Mean 
Width of Primary Shear Deformation Zone. 
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4,6. Eff.C of cutting V.locity on ‘f;* 

M.on Width o( Primary Sh.or D.tormotion 
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CHAPTER 

\ 

c 

CONCLUSION AtID SUGGESTIONS FOR FUTURE WORK 


CONCLUSION t 

From the results obtained, in the present work, the 

following conclusions can he drawn* 

(i) Using the^ompressiblej non»Newtonian fluid flow for 
describing the plastic deformation during metal cutting,f 
the primary shear deformation a one can he calculated 
reasonably accxirately, 

(ii) A cut off value of 35^* of the maximum shear strain rate 
invariant can he used for approximately demarcating 
the primaiy shear zone, 

(ill) PEM can be used for solving Idae velocity and strain 
rate distribution In the primary shear a on® during 
cutting. It is recommended that Penalty Formulation 
used (which ellirainates the pressure variabl®) for 
solving the flow equation* Since inclusion of pressure 
leads to serious numerical difficulty, 

3TOE3TI0NS FOR FUTURE WORK t 

1 , The effects of temperature on plastic flow properties in 
energy equation and flow equation must he considered. 

The analysis can he generalised for oblique cutting usini 
three -dimensional formulation, 

3t Different type of yiacosity equations can he tried for 
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